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The properties of exotic nuclei are the focus of the present research. Two-neutron halo structures of neutron- 


rich '’:!°B were experimentally confirmed. We studied the formation mechanism of halo phenomena in 


17,198 


using the complex momentum representation method applied to deformation and continuum coupling. By 
examining the evolution of the weakly bound and resonant levels near the Fermi surface, s—d orbital reversals 
and certain prolate deformations were observed. In addition, by analyzing the evolution of the occupation 
probabilities and density distributions occupied by valence neutrons, we found that the ground state of !°B did 
not exhibit a halo and the ground states of 17B and !°B exhibited halos at 0.6 < 82 < 0.7 and 0.3 < B2 < 0.7, 
respectively. The low-! components in the valence levels that are weakly bound or embedded in the continuous 


spectrum lead to halo formation. 
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I. INTRODUCTION 


Exotic phenomena, including halos, have attracted consid- 
erable interest from physicists [1—4] and are one of the most 
popular frontier topics in radioactive nuclear beam physics. 
In isotope chains, with an increase in the neutron number, 
nuclei near the neutron-rich dripline may remarkably appear 
as neutron halos. As ‘Li was discovered to be an s-wave 
two-neutron (2n) halo nucleus [5], the formation of a halo 
is considered to be a large occupation of low-orbital-angular- 
momentum (l = 0, 1) levels close to the threshold of particle 
emission in weakly bound systems [6—8]. Research on halo 
structures, both theoretically and experimentally, is of great 
importance in understanding the structure of the nucleus. 

On the experimental side, extensive research on halo nuclei 
with small mass numbers has led to a clear understanding of 
their halo structures; for example, 8B is a proton (1p) halo nu- 
cleus [9]; ĉHe [10, 11], ‘Li [5, 12], and “Be [13, 14] are 2n 
halo nuclei; and HBe [15, 16] is a neutron (1n) halo nucleus. 
Subsequently, ??C, !”B, and 19B were found to be s-wave 2n 
halo nuclei with large matter radii [17, 18], a small percent- 
age of s-orbital spectroscopic factors [19], and an enhanced 
electric dipole strength [20], respectively. In the heavier mass 
region, 37Mg [21] (the heaviest halo nucleus observed thus 
far) and 3!Ne [22] and ?9F [23] (the heaviest Borromean nu- 
cleus to date) were found to be p-wave halos at or near the 
inversion island by single-neutron (17) removal reactions and 
reaction-cross-section measurements, respectively. 

On the theoretical side, the search for new 2n-halo sys- 
tems [24] and new halo features, such as core—halo shape de- 
coupling [25, 26] and the Efimov state [27, 28], have also 
attracted considerable interest. It can be seen that valence 
neutrons predominantly occupy p and s orbitals, so that cor- 
relations between f—p or s—d orbitals play an essential role in 
the formation of halo phenomena. These strong quadrupole 
correlations between f—p or s—d orbitals can drive the sys- 
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tem to break spherical symmetry and undergo deformation, 
allowing the number of single-particle levels with low-l com- 
ponents to increase, resulting in abundant halos at the neu- 
tron dripline [29]. Specifically, Hove et al. [30] used the 
example of a medium-heavy nucleus ’*Ca and found that the 
presence of a slight, long-range s-wave tail is also sufficient to 
bring about halos in very weakly bound neutron-rich systems. 
Therefore, it is necessary to study the effects of quadrupole 
deformation and continuum effects on the halo structure of 
exotic nuclei. 

In this study we investigated the neutron halo structures of 
B isotopes. 15B and 17B neutron halo structures have been 
reported [31, 32]. For !°B, the 2n separation energies (S2n) 
were 3.747 MeV [33], the quadruple deformation length 6 
was 1.17 [34], the quadruple deformation parameter G2 was 
0.61 [35], and the percentage of the s orbital was 63% [36]. 

For 17B, a strong mixing of s and d waves was observed 
based on the energy dependence observed by the interaction 
cross section [37]. Tanihata [24] considered large amplitudes 
of 28, /2 waves for all B isotopes with neutron number n > 8. 
For 17B, Sən is 1.37 MeV [38], 6 is 1.2 [34], Bo is ~0.6 [39] 
or 0.56 [35], the neutron skin thickness is 0.51 (0.11) fm [40], 
and the neutron has a large matter radius [40, 41]. The per- 
centages of s orbital in 17B are 36(19)% [31], 69(20)% [32], 
50(10)% [37], and 53(21)% [42], respectively. In 2021, Yang 
et al. [19] obtained spectral factors for the 1s1/2 and 0d5/2 
orbitals, with the 1s,/2 orbitals contributing only 9(2)% of 
the spectral factor, which is much smaller than the values 
given in Refs. [31, 42], revealing that 17B has a well-defined, 
but not dominant, neutron halo. The structure of this weak 
halo has been explained by the deformed relativistic Hartree— 
Bogoliubov theory in the continuum (DRHBc) [19], indicat- 
ing the significance of deformation effects. 

19B is the heaviest bound nucleus of the B isotopes, and, 
because 18B is unbound, 19B is the Borromean nucleus, sug- 
gesting a 2n halo structure. 19B has a particularly low two- 
neutron separation energy (S2n) of 0.089" 9:560 MeV [38], 
enhanced interaction cross section [31], and large deforma- 
tion [35, 43, 44]. Direct mass measurements revealed that the 
valence neutrons of !9B are dominated by d waves, which 


suppress halo formation [31, 42, 45]. However, Cook et 
al. [20] found a significant 2n halo in 1°B using Coulomb 
dissociation. Recently, !°B has been suggested as a deformed 
2n halo candidate with shape-decoupling effects [35]. The 
structures of !°B are also related to the newly discovered un- 
bound isotopes 7°B and ?!B. 

However, the deformation of B isotopes close to the neu- 
tron dripline has not been experimentally determined. The 
shape of the nucleus is determined by the shell structure of 
protons and neutrons. Therefore, it is of interest to study 
quadrupole deformation and continuum effects on halo for- 
mation in B isotopes. 

The continuum and particularly the resonant states in the 
continuum play a crucial role in the formation of exotic 
phenomena. Therefore, it is essential to reasonably handle 
the resonant states in the continuum, especially those near 
the threshold. Given the significance of resonant states in 
studying exotic phenomena, physicists have developed many 
theoretical resonance models, such as the Green’s function 
method [46], the complex scaling method (CSM) [47], and 
the complex momentum representation (CMR) method [48]. 

To study the halo phenomenon, we used the CMR method 
based on the coupled-channel approach to solve the equations 
of motion of the deformed nucleus. It is worth noting that 
the quadrupole—quadrupole force in the spherical mean field 
is included in the deformation mean field and would make 
the single-particle figure work much better in the deforma- 
tion nucleus than in the spherical nucleus. To determine the 
appearance of the halo phenomenon more accurately, it is es- 
sential to use models in which the radial wavefunctions of the 
weakly bound particles are treated appropriately in this study. 

The CMR method can be used to obtain the bound states, 
resonant states, and continuum on an equal basis [48, 49]. 
This method is not only very effective for narrow resonances 
but also can be reliably applied to broad resonances. In par- 
ticular, the results of the calculations have been proven to 
be independent of any nonphysical parameters, which en- 
sures that high-precision energies and widths can be obtained. 
The CMR method is widely used to study the properties of 
weakly bound nuclei [50, 51]. We have recently studied 
in detail the structures of exotic spherical nuclei using the 
CMR method combined with the relativistic meson exchange 
model [48, 52] and relativistic point coupling model [53]. 
Furthermore, pairing correlations are treated with the BCS 
approximation, coupling bound states to resonant states and 
describing exotic phenomena such as halos and giant halos 
very well [54, 55]. 

In addition, we extended the CMR method to a deformed 
system. We have systematically investigated experimen- 
tally discovered and undiscovered halo structures, establish- 
ing that *4Ne [50], 37Mg [49], 53Ar [56], 34Na [57], and 
29,.31F (58, 59] (or s waves or a mixture of s and p waves) 
are p-wave halo nuclei and that °C [60], "5Cr [56], and 
TTFe [56] are s-wave halo nuclei. Recently, Wang et al. in- 
vestigated neutron-rich ’°Cr using the CMR method in a rela- 
tivistic mean-field framework, which applies to deformed nu- 
clei. They found that Cr is a d-wave deformed halo nu- 
cleus [61]. The ground-state properties of Ti isotopes were 


also investigated using the relativistic mean-field framework 
CMR method [62]. The predictions of halo phenomena in the 
neutron-rich region close to the neutron dripline will help to 
experimentally explore the deformation halos of more nuclei 
in medium-mass regions. The present research demonstrates 
the advantages of the CMR method, which will help us ex- 
plore further deformed halo nuclei using the CMR method. 

In this study, we systematically investigated the exotic phe- 
nomena of 15:17-19B using the CMR method. The energies 
and widths of the single-particle bound and resonant states, 
component ratios of the orbitals occupied by valence neu- 
trons, and radial density distribution were obtained. The re- 
sults revealed that '”!9B may be prolate halo nuclei, which 
can be well understood because of the neutron halo phe- 
nomenon. In other words, observing the halo phenomena in 
17,19B can be used as critical evidence for the prolate defor- 
mation of '”+9B ground states. The theoretical formalism is 
described in Sect. II. Numerical details and results are pre- 
sented in Sect. II. Finally, a summary is provided in Sect. 
IV. 


Il. FORMALISM 


We briefly introduce the theoretical framework of the CMR 
method to explore the ground-state properties of neutron-rich 
nuclei. The details of the CMR method can be found in 
Refs. [48, 50]. The main formulas of the CMR method are 
provided for the convenience of later discussion. We intro- 
duce a theoretical framework using the Hamiltonian 


H=T +V, (1) 


where T = p? /2m denotes the kinetic operator with momen- 


tum p = hk (where k is the wavevector). Following Ref. [63], 
the interaction potential V consists of the central potential, 
deformation potential, and spin-orbit coupling potential: 


Veent (r) = V (r) , 
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where A is the reduced Compton wavelength of the nucleon, 
h/m,c, and 
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The diffuseness a = 0.67 fm, the radius R = rọ A?/? with 
ro = 1.27 fm, and the depth of potential Vọ = —51+334>4 
MeV depends on the mass number A, proton number Z, and 
neutron number N. To obtain the bound states, resonant 
states, and continuum on an equal footing, the Schrödinger 
equation is transformed from a coordinate representation to a 
momentum representation as 


J dË (KIHE yE) = EY), 8) 


where w(k) are the momentum wavefunctions. For axially 
deformed nuclei, the third component of the total angular mo- 
mentum Q and parity m are suitable quantum numbers. The 
wavefunctions o (K) can be expanded as 


valk = a ss 


where f'i (k) and ¢1;9(Qx) are the radial and angular com- 
ponents, respectively. By substituting Eq. (4) into Eq. (3), 
the Schrödinger equation becomes 
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The above Schrödinger equation is solved in complex mo- 
mentum space by diagonalizing the Hamiltonian matrix. For 
the sake of computational simplicity, the integral of Eq. (5) is 
transformed into summation form, and Eq. (5) becomes the 
matrix equation. In this study, it was combined with the cou- 
pled channel method to describe the deformed nuclei. Details 
can be found in Refs. [59, 60]. 
For every Nilsson level, the wavefunction is extended with 
a set of spherical components. The wavefunctions in the co- 
ordinate representation can be obtained using the following 
transformation: 
vo (r) = ——_ J ke val) 
(27)? /2 


By assuming that Yo (r) holds the form 


be (F = LO) (r) b1j0(Qr), (6) 
where 
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the corresponding radial density distributions can be calcu- 
lated as 
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Fig. 1. (Color online) RMS radius of single-particle levels occupied 
by the last valence neutron as a function of the mass number A for 
13-19B | The corresponding deformation $2 and the quantum num- 
bers of the level occupied by the last valence neutron are displayed 
in every subfigure. 


Til. NUMERICAL DETAILS AND RESULTS 


The details of the numerical calculations are the same 
as those in Refs. [59, 60], where the number of coupling 
channels, N,, momentum truncation maximum value kmaz» 
and number of Gauss—Legendre grids, N; are taken to be 
8, 4.0 fm™t, and 120, respectively. The CMR calculations 
achieved the required precision under these numerical con- 
ditions. Because quadrupole deformation and continuum ef- 
fects play an essential role in the formation of the halo of 
exotic nuclei [35, 58, 64], the quadrupole deformation param- 
eter G2 was taken from —0.4 to 0.7 in steps of 0.02. Based 
on Refs. [48, 50], we found that the positions of the reso- 
nant states do not depend on the integration path. There- 
fore, we chose a sufficiently large contour (kı = 0.0 fmt, 
ko = 0.5— 10.5 fm7!, kg = 1.0 fmt, and kmar = 4.0 fm7 1) 
to expose all relevant resonances. Based on the numerical de- 
tails presented, we first reveal the contributions of quadrupole 
deformations and continuum effects to the ground-state prop- 
erties of 17B. 

The nuclear radius is a significant physical quantity that 
describes the exotic properties of nuclei. It is known that the 
root mean square (RMS) radius of halo nuclei (nuclei that 
form a halo away from the periphery of the core) increases 
suddenly compared to that of the neighboring nucleus. It can 
be seen from the definition of halo nuclei that the RMS radii 


of the orbits occupied by the last few valence neutrons are 
of crucial importance in researching halo phenomena. There- 
fore, we calculated the RMS radii for all the single-particle 
levels in the considered deformation range. For convenience, 
as shown in Fig. 1, we show only the RMS radii of the levels 
occupied by the last two valence neutrons of !*-!°B for the 
quadrupole deformations G2 = —0.1, 0.0, 0.1, 0.2, 0.3, 0.4, 
0.5, 0.6, and 0.7. 

The halo structure was experimentally confirmed for 
17,19B [19, 20]. The last two valence neutrons of !”B oc- 
cupy the levels 5/2[202], lds /2, 1/2[220], and 3/2[211] at 
Bo = —0.1, Bo = 0.0, Bo = 0.1-0.2, and Bo = 0.3- 
0.7, respectively (see Fig. 2). We can see that the RMS ra- 
dius of level 3/2[211] occupied by the last two valence neu- 
trons of !”B significantly increases compared to that of '°B at 
B2 = 0.6 and 2 = 0.7. However, we noticed that the RMS 
radius of ‘”B increases slightly at 82 = 0.1-0.2, as a result of 
the valence neutrons occupying the resonant level 1/2[220]. 
This deformation region G2 = 0.1-0.2 was excluded because 
17B is a weakly bound nucleus. Similar phenomena were ob- 
served in !°B, where the RMS radius increased significantly 
compared to that of '7B at 8. = 0.3-0.7. For !°B, there 
was no significant increase in the RMS radius over the range 
of deformations compared to the neighboring nuclei. Simi- 
lar to the phenomena at 62 = 0.1-0.2 in 17B, the increased 
RMS radius for 15B at 82 = —0.1-0.2 is caused by the last 
two valence neutrons occupying an extremely weakly bound 
level 1/2[200]. However, the deformation region G2 = —0.1- 
0.2 was excluded because the two-neutron separation energy 
of the weakly bound level 1/2[200] was much less than the 
experimental value of 3.747 MeV [33]. These phenomena 
suggest that '7B and !°B may be prolate halo nuclei with de- 
formations 82 = 0.6-0.7 and 82 = 0.3-0.7, respectively. In 
comparison, !°B is only a prolate deformation nucleus, which 
is relatively close to the results obtained from DRHBc calcu- 
lations [35]. 

The neutron halo structure can be identified by analyzing 
the valence neutron levels. In Fig. 2, the Nilsson diagram 
of 17B, including the resonant and bound levels, is shown 
to better study the effects of deformation and the contin- 
uum on halo structures. The bound and resonant states are 
marked with solid and dashed lines, respectively, and the cor- 
responding Nilsson orbitals are marked with quantum num- 
bers Q[/Nn-A]. Levels 1p: /2 and 251/2, and 1d5/2 and 1d3/2 
have large shell gaps between them when spherical; that is, in 
the spherical case, the old magic number N = 8 is verified, 
and the new magic number N = 16 is suggested. 

For the spherical configuration (82 = 0), level 2s1/2 is 
lower than lds /25 which means s—d orbitals are inverted com- 
pared to the level order in the stable nuclei. This is similar to 
the case of tB from DRHBc calculations [35]. In the spher- 
ical case, the resonant orbital 1d5 2 is occupied by the last 
two valence neutrons of !7B, which have a large centrifugal 
barrier and negative two-neutron separation energy, suppress- 
ing the formation of halos. That is, the spherical shape does 
not support the formation of the halo phenomenon. Because 
the halo phenomena of !”B have been observed experimen- 
tally [19], this suggests that 17B is not a spherical nucleus but 
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Fig. 2. (Color online) Neutron single-particle levels of "B as a func- 
tion of G2. Every Nilsson orbital is labeled with the quantum num- 
bers Q[Nn-A]. The bound and resonant levels are marked by solid 
and dashed lines, respectively. 


rather a deformed halo nucleus. Whether 17B is a prolate or 
oblate deformation halo nucleus is an problem to be solved. 


Driven by the deformation, the spherical levels split and 
the degenerate single-particle levels become undegenerated. 
We can see that some of these levels split off, some fall 
with deformation, and some rise with deformation. The reso- 
nant state 1ds5/2 splits into Nilsson levels 1/2[220], 3/2[211], 
and 5/2[202], and the weakly bound level 2s, /2 becomes 
1/2[200]. 

Nilsson levels 1/2[200], 5/2[202], 1/2[220], and 3/2[211] 
are occupied by the last two valence neutrons of !’B at 
—0.4 < Bo < —0.14, —0.14 < B2 < 0,0 < B2 < 0.26, 
and 0.26 < B2 < 0.7, respectively. Given that 17B is 
a weakly bound nucleus, the last valence neutron may oc- 
cupy the weakly bound level 1/2[200] with a deformation 
—0.4 < By < —0.14 or the weakly bound level 3/2[211] 
with a deformation 0.26 < 82 < 0.7. 


On the oblate side, when pairing correlations are not con- 
sidered, Nilsson level 1/2[200] is occupied by the last two 
valence neutrons of 17B at —0.4 < Ba < —0.14. It should 
be indicated that the pairing correlations play important role 
at halo formation [65, 66]. Given the impact of pairing cor- 
relations, the last two neutrons can partially occupy level 
1/2[200]. The two-neutron separation energies of the weakly 
bound level 1/2[200] are <0.75 MeV, which is much lower 
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Fig. 3. (Color online) Same as Fig. 2, but for the widths for all the 
resonances of interest. 


than the experimentally obtained value of 1.37 MeV [38]. 
This phenomenon suggests that 17B is not an oblate nucleus. 


In contrast, when pairing correlations are not considered, 
the Nilsson level 3/2[211] is occupied by the last two va- 
lence neutrons of !"B at 0.26 < B2 < 0.7. The Sən value 
of level 3/2[211] reaches 1.3-1.6 MeV with a deformation 
0.6 < B2 < 0.7, which is close to the experimentally ob- 
tained 1.37 MeV [38] and, when combined with the analysis 
in Fig. 1, suggests that favorable conditions for halo forma- 
tion are met in this deformation region. Given the impact of 
pairing correlations, the last two neutrons can partially oc- 
cupy the level 3/2[211]. This implies that !"B is a prolate 
deformation halo nucleus. This does not contradict the pre- 
diction of a prolate deformed halo nucleus for !B. 


Figure 3 shows the widths of the single-particle resonant 
states in 17B. For simplicity, the widths in Fig. 3 are the same 
as those for the levels shown in Fig. 2. Similar to the levels, 
there is a shell structure in the widths. Over the deformation 
range under consideration, the gap at N = 16 is relatively 
large, not only in a spherical shape but also in a deformed 
shape, which means that a new magic number N = 16 is 
supported in !"B. The width of the resonant state 3/2[211] 
is consistently small and lies near zero potential energy. As 
62 increases, the width of resonant state 3/2[211] shrinks and 
the state eventually disappears, which is due to the resonance 
state 3/2[211] becoming a bound state. Therefore, it is rela- 
tively easy to form a halo at the level 3/2[211] in 17B. The 
experimentally detected lifetime of !7B will help us better 
understand the halo structure of 17B because the widths are 
inversely proportional to the decay lifetime. 

The above discussions indicate that the last two valence 
neutrons of 7B occupy the Nilsson level 3/2[211] and that 
17B is a prolate halo nucleus with a deformation 0.6 < 62 < 
0.7. To reveal the spherical component of the Nilsson orbitals 
occupied by the last valence neutron, the occupation probabil- 
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Fig. 4. (Color online) Occupation probabilities of major configu- 
rations as a function of 82 for the single-particle level 3/2[211] in 
17 
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Fig. 5. (Color online) Wavefunctions of every component for the 
single neutron orbital 3/2[211] at 82 = 0.64 for "B. 


ities of the level 3/2[211] are shown in Fig. 4. It can be seen 
that the configuration d5/2 exceeds 90%. Configuration d3/2 
increases slowly, reaching ~6% at 32 > 0.62. This suggests 
that the Nilsson level 3/2[211] is composed mainly of the d 
configuration, which is consistent with the results obtained in 
Ref. [67]. Meanwhile, the occupation probabilities of level 
1/2[200] were plotted (and, as these are similar to those in 
Fig. 7, they are not shown here). The occupations of config- 
urations 51/2 and d5/2 remain almost constant at 62 > 0.3, 
both reaching ~47%, indicating that level 1/2[200] is mainly 
composed of configurations ds /2 and s1/2. This suggests that 
the occupation of level 1/2[200] makes an important contri- 
bution to the distribution of diffuse density in 17B. 


To better understand the effects of quadrupole deformation 
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Fig. 6. (Color online) Radial density distributions multiplied by r? 
for the single-particle orbitals 3/2[211], 1/2[200], and 1/2[101] with 
b2 = 0.64 in "B. 


and single-particle levels occupied by valence neutrons on the 
halo structure, as shown in Fig. 5, we plotted the radial wave- 
function for the 3/2[211] orbital at G2 = 0.64. As in Fig. 4, 
the radial wavefunctions are shown only for the real part of 
the configurations d3/2, ds 2, 97/2, and gg/2. The other spher- 
ical components and the imaginary parts can be ignored. The 
dependence of the wavefunctions of configurations d3/2 and 
ds 2 (l = 2) on the radius r at long distances is not signif- 
icantly different from that in the bound state. However, the 
wavefunctions of configurations d3/2 and d5/2 have large ra- 
dial density distributions at r ~ 3 fm. It is well known that 
the configuration d-wave configuration does not support the 
formation of halos owing to its large centrifugal potential bar- 
rier. Therefore, whether the halo phenomenon in 17B occurs 
at 0.6 < 82 < 0.7 requires further exploration. 

The exotic characteristic of the halo is its density distri- 
bution with a long tail compared to its neighboring levels. To 
check whether the levels 3/2[211] are occupied by the last two 
valence neutrons of !7B at 0.6 < B2 < 0.7 is a reasonable ap- 
proach. As shown in Fig. 6, we explored the density distribu- 
tion of the single-particle orbital 3/2[211] at 82 = 0.64. For 
comparison, the radial density distributions of levels 1/2[200] 
and 1/2[101] are also shown. The radial density distributions 
of the bound levels 1/2[200] and 1/2[101] did not diffuse at 
long distances. However, the weakly bound level 3/2[211] 
can be found trailing a longer tail at greater distances. Ac- 
cording to Fig. 2, we can see that the level 3/2[211] decreases 
slowly at 0.6 < 62 < 0.7. These phenomena indicate that 
the last two valence neutrons of !”B occupy level 3/2[211] 
when 0.6 < B < 0.7 and weak halo phenomena are ob- 
served. These results are consistent with the halo observed for 
17B [19]. In summary, the observed halo phenomena suggest 
that 7B is a deformation halo nucleus with a deformation 
parameter 0.6 < 8 < 0.7 and that its ground state consists 
mainly of the d configuration. However, we did not obtain 
the exact percentage of the 2s, /2 orbital, which also shows 
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Fig. 7. (Color online) Same as Fig. 4, but for 15B with the single- 
particle level 1/2[200]. 


that the prerequisite for the occurrence of neutron halos is 
not the dominance of components with smaller orbital angu- 
lar momentum. 

The reliability of the CMR method for studying exotic phe- 
nomena is well demonstrated by investigating the halo struc- 
ture in 17B. Hence, to obtain the ground-state properties of 
the other neutron-rich B isotopes, we further explored the ex- 
otic structures of neutron-rich 15B and !°B using the CMR 
method. 

Because !°B has only two fewer neutrons than !’B, the or- 
der of the single-particle levels does not change, except for a 
slight overall upward shift. When pairing correlations are not 
considered, as shown in Fig. 2, it can be seen that the last two 
valence neutrons of !°B occupy levels 5/2[202] and 1/2[200], 
respectively. The Sə» value of the level 5/2[202] is <1.56 
MeV at 62 > —0.4, and the Sən value of level 1/2[200] is 
<3.4 MeV and >4 MeV at (2 > 0.5 and 82 < 0.56, respec- 
tively, which is not in agreement with the experimentally ob- 
tained 3.747 MeV [33]. These phenomena suggest that 15B 
may be a deformed nucleus, with 82 between 0.5 and 0.56. 
Given the impact of pairing correlations, the last two neutrons 
can partially occupy level 1/2[200]. This does not contradict 
the prediction of a prolate deformed nucleus for '°B. 

To confirm the results obtained above, the occupation prob- 
abilities of levels 1/2[200] are plotted in Fig. 7. The oc- 
cupations of configurations 51/2 and ds5/2 decrease and in- 
crease rapidly with an increase in 82, respectively. When 
B2 > 0.3, their occupation remains almost constant, both 
reaching ~50%, indicating that the level 1/2[200] is mainly 
composed of configurations s1/2 and ds 2. 

As in Fig. 5, we show the radial wavefunction at 62 = 0.52 
for the single-particle orbital 1/2[200] of 15B in Fig. 8. The 
wavefunction of configuration 51/2 decreases gradually as r 
increases from 0 to 3 fm. The wavefunction of configuration 
ds/2 has a large radial density distribution at r ~ 3 fm. As 
r increases further, the wavefunctions of configurations sj 2 
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Fig. 8. (Color online) Same as Fig. 5, but for B with the single- 
particle level 1/2[200]. 
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Fig. 9. (Color online) Same as Fig. 6, but for **B at 82 = 0.52. 


and ds/2 converge rapidly to zero at long distances. 

Meanwhile, in Fig. 9, we plotted the density distributions 
of level 1/2[200] in 1°B at Bg = 0.52. For comparison, 
we also show the density distribution for the 3/2[211] and 
1/2[101] levels. The density distribution of the bound level 
1/2[200] converges rapidly at longer distances, exhibiting no 
halo tails, suggesting that halo formation is not supported 
when the valence neutrons occupy level 1/2[200]. That is, 
15B is a prolate nucleus, but no halo phenomenon exists. 

We performed a similar investigation on 1°B. Although 
19B has only two more neutrons than 17B, the contribution 
of single-particle levels near the continuum threshold to the 
exotic phenomena becomes more significant, in addition to 
a slight downward shift at the single-particle level. This is 
consistent with the results obtained in Refs. [35](in Fig. 2). 
When pairing correlations are not considered, the Nilsson 
levels 3/2[211], lds /2, 3/2[211], and 1/2[220] are occupied 
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Fig. 10. (Color online) Same as Fig. 6, but for '°B at G2 = 0.36. 


by the last two valence neutrons of !9B. Given that !°B is a 
weakly bound nucleus, the two-neutron separation energy is 
0.0897} 3g, MeV [38]. Therefore, at 0.3 < 8 < 0.7, level 
1/2[220] can be occupied by the last two valence neutrons 
of '°B. Given the impact of pairing correlations, the last two 
neutrons can partially occupy level 1/2[220]. This is not con- 
tradictory to the prediction of a prolate deformed halo nucleus 
for !°B. However, whether halo phenomena appear in this de- 
formation region requires further exploration. 

To verify whether quadrupole deformation with 0.3 < 
B < 0.7 is a reasonable conjecture in 1B, we explore the 
density distribution of the single-particle orbital 1/2[220] at 
82 = 0.36 in Fig. 10. For comparison, the radial density 
distributions of levels 5/2[202], 3/2[211], and 1/2[200] are 
also shown. The radial density distributions of the bound lev- 
els 3/2[211] and 1/2[200] were not very diffuse at large dis- 
tances. It is noteworthy that the level 1/2[220] has a larger 
diffuse density distribution than the resonant level 5/2[202]. 
These phenomena indicate that the last two valence neutrons 
of !°B occupy level 1/2[220] when 0.3 < 8 < 0.7 and halo 
phenomena are observed. These results are consistent with 
the observed halo [20] and the predicted deformation [35]. 

In addition, we also provide the occupation probabilities of 
level 1/2[220] in Fig. 11. This level consists mainly of the 
configurations 1/2, d3/2, and ds/2, whereas the other con- 
figurations are negligible. Configuration s,/2 is dominant 
compared to configurations d3/2 and d5/2. The occupation 
probability of configuration 51/2 increases rapidly and then 
remains almost constant as (> increases; it still exceeds 65% 
on the prolate side; that is, the wavefunction of the weakly 
bound level 1/2[220] is composed mainly of configuration 
81/2, which supports the formation of the s-wave halo in '°B. 

To better understand the effect of single-particle levels 
with low orbital angular momentum on the halo structure, in 
Fig. 12 we plotted the radial wavefunction for orbital 1/2[220] 
at 62 = 0.36. The radial wavefunctions are shown here only 
for the real and imaginary parts of configurations 51/2, d3/2, 
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Fig. 11. (Color online) Same as Fig. 4, but for ‘°B with the single- 
particle level 1/2[220]. 
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Fig. 12. (Color online) Same as Fig. 5, but for ‘°B with the single- 
particle level 1/2[220]. 


and ds 2. The wavefunctions of the real and imaginary parts 
of the configuration s;/2 (l = 0) extend to a considerable 
range in the coordinate space and exhibit significant halo phe- 
nomena. In contrast, the dependence of the wavefunctions of 
the real and imaginary parts of the configuration d5/2 (l = 2) 
on radius r is not significantly different from that in the bound 
state. Although the wavefunctions of the real and imaginary 
parts of configuration d;/2 have a large radial density distri- 
bution at r ~ 3 fm, these local wavefunctions do not support 
halo formation because of the large centrifugal potential bar- 
rier and small occupation probabilities. The halo phenomena 
observed in '°B suggest that '°B is a deformation halo nu- 
cleus with a deformation parameter 0.3 < 6 < 0.7 and that it 
is the ground state consisting mainly of the s configuration. 
The present method allows us to identify 15B as a possible 
prolate nucleus with 32 between 0.50 and 0.56, and the halo 


phenomena in +B and !°B may occur at quadrupole defor- 
mations 0.6 < 2 < 0.7 and 0.3 < 6 < 0.7, respectively. 
The above results indicate that the presence of !7B and !°B 
halo structures can be qualitatively explained in a valid model 
that considers quadrupole deformation and continuum spec- 
tral effects, despite its various limitations. The ground-state 
structures of 15B are mainly the result of a certain degree of 
deformation. However, the ground-state structures of 17B and 
19B are primarily the result of the interaction between s and 
d orbitals, large quadrupole deformation, and continuum ef- 
fects. 


IV. SUMMARY 


Inspired by the recently reported experimental work on the 
halo phenomena in '7:!°B [19, 20], the role of quadrupole 
deformation and continuum effects on the halo structure of 
17B was studied using the CMR method. 

Single-particle orbitals, including the bound and resonant 
states, were obtained, and an inversion of the s—d orbitals was 
found to be a critical factor in yielding the halo phenomena. 
At 2 = 0, the neutron orbital 281 /2 is lower than that of 
1d5/2, and the valence neutrons are mainly dominated by the 
d wave. Therefore, !"B is not a spherical nucleus. Driven 
by the deformation, the ground state of !’B is prolate (0.6 < 
Ba < 0.7) and a halo appears. The calculated occupation 
probabilities, radial density distributions, and wavefunctions 
support the results obtained by using the Nilsson levels. This 
indicates that the ground state of 17B is a prolate halo nucleus, 
with 0.6 < b2 < 0.7. 

As shown in Figs. 7-12, the ground state of '°B also ex- 
hibits a large deformation but no halo phenomenon. By ob- 
serving the radial density distribution and the configurations 
of the orbitals occupied by the valence neutrons, we find that 
the ground state of !°B appears as halos at 0.3 < By < 0.7 
and that the ground state consists mainly of the s configura- 
tion. 

In this study, we used the CMR method to present some ba- 
sic physics in a simple manner. By studying the quadrupole 
deformation and continuum effects in the ground-state prop- 
erties of 15B, 17B, and !°B, we found that the ground-state 
structure in !°B may be significantly influenced by defor- 
mation effects and that the halo structures in !7-!°B may re- 
sult from the interaction between s—d orbital reversals, large 
quadrupole deformations, and continuum effects. This work 
shows that further exploration of the relationship among halo 
features, deformation effects, and coupling of continuum 
spectra in more exotic nuclei is essential and interesting. Fur- 
thermore, the development of a CMR model is important, and 
future work will consider pairing correlations in this theoret- 
ical formalism. 
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